Electrically controlled persistent spin currents at the interface of multiferroic oxides 
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We predict the appearance of a persistent spin current in a two-dimensional electron gas formed 
at the interface of multiferroic oxides with a transverse helical magnetic order. No charge current is 
generated. This is the result of an effective spin-orbit coupling generated by the topology of the oxide 
local magnetic moments. The effective coupling and the generated spin current depend linearly on 
the magnetic spiral helicity which, due to the magneto-electric coupling, is tunable by a transverse 
electric field offering thus a new mean for the generation of electrically controlled persistent spin 
currents. 



INTRODUCTION 



II. THEORETICAL FORMULATION 



Nanoscience research is fueled by the spectacular func- 
tionalities that emerge from the controlled composition 
of different materials down to the atomic level. A 
recent example is the appearance of a metallic phase 
with a high carrier mobility confined to the inter- 
face between insulating oxides 1 such as LaTiOs/SrTi03 
or LaAlOa/SrTiOa 2 . This sheet of two-dimensional 
electron gas (2DEG) has been laterally confined and 
patterned 3 to achieve nanometer-sized tunnel junctions 
and field-effect transistors 4 , thus paving the way for 
oxide-based Nanoelectronics 3 with a multitude of tech- 
nological applications 5 . Further functionalities are ex- 
pected when utilizing the residual properties of the ox- 
ides. E.g., an important group of Mott insulating oxides 
such as RMn0 3 (R= Tb, Dy, Gd, and Eui-^) 6 and 
LiC^O^ are multiferroics with a noncollinear magnetic 
phase. The origin of the spontaneous electric polariza- 
tion is argued 8 to be the spin current associated with 
the spiral magnetic order. As shown experimentally, due 
to the magneto-electric coupling, the helicity associated 
with the spin spiral structure of the multiferroics is tun- 
able from clockwise to counterclockwise type by a small 
electric field (~ lkV/cm) 9 . 

In this paper we show theoretically, that a 2DEG 
formed at the surface of a multiferroic oxide (FigJJ) such 
as the ab plane of TbMn03 6 experiences an effective spin- 
orbit interaction (SOI) that linearly depends on the car- 
riers wave vector and on the helicity of the oxide's mag- 
netic order and hence is controllable by a lateral electric 
field. As a result an electrically tunable persistent spin 
current is shown to build up in the 2DEG. No charge 
current is generated. The origin of this effect lies in the 
topological structure of the local magnetic moment at the 
oxides interface. Spin currents are actively discussed in 
the field of semiconductor-based spintronic a 10 ! 11 ! 12 ! 13 ! 14 . 
There, SOI plays also a vital role. In semiconductors 
however, a finite dissipative charge current is also gen- 
erated by the applied in-plane electric field. Hence, the 
persistent spin current in insulator 15 , as uncovered here, 
has a decisive advantage, as compared with metals^ and 
semiconductors 12 and adds a new twist to oxide electron- 
ics. 



As sketched in FigHJ we consider a 2DEG, as realized 
in Refs.[1.2.3.4] at the interfaces of oxides layers, how- 
ever one of the layer should be a spiral multiferroic oxide 
such as TbMn03 6 or LiC^O^. The spiral structure de- 
fines the x — z plane, whereas the 2DEG is confined to 
the x — y plane (cf. Fig{T]). At low temperature, the ox- 
ide local spin dynamics is much slower than the 2DEG 
carrier dynamics and hence we can treat the oxide local 
moments as classical and static. A carrier in the 2DEG 
with a charge e experiences an effective (real) internal 
magnetic field due to the magnetic spiral and the em- 
bedding mediu m 17 i 18 i 19 which results in a nonlocal vector 
potential A in . The effect of A in on the charge carriers 
dynamics is subsidiary compared to that of the exchange 
field Jn r where n r is a local unit vector field describing 
the geometry of the localized magnetic moments at the 
oxides interface and J is the coupling strength. The ex- 




FIG. 1: (Color online) Schematic of proposed system. The 
£-axis is defined along the direction of spiral ordering, for 
example, the [110] direction in TbMnOs. The spiral plane of 
multiferroic (below) is perpendicular to the 2DEG (above). 
Due to magneto-electric coupling the spin helicity is controlled 
by the transverse electric field E z . 
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change interaction originates from the Coulomb repulsion 
on the localized moments and from Hund's-rule coupling 
in the magnetically ordered phase & 2Q i 21 . Thus, the sin- 
gle particle dynamics in the 2DEG is governed by the 
Hamiltonian 5 



H = h k 



hj 



2rn 



Jn r ■ a 



(i) 



where m is the effective electron mass, a is the vector 
of the Pauli matrices and P is the momentum operator. 
n r is given by the local magnetization at the multiferroic 
surface, i.e. n r = (sin# r , 0, cos# r ) where r = q m -r with 
q m = (g,0,0) being the spin-wave vector of the spiral. 
A in is not included in Eq.(pp).— 

Applying the unitary local gauge transformation in 
the spin space U g — exp(—i6 r <j y /2), the spatially non- 
homogeneous term hj is transformed into the diago- 
nal term 24 hj = U^hjU g = Ja z (Hereafter, trans- 
formed quantities are marked by a tilde). Physically, 
this amounts to a rotation of the local quantization 
axis to align with n r at each site. a y = a y because 
[U g , <J y ] =0. We find further d x = (<j x cosO r — cr z sin# r ) 
and d z = (a x sm6 r + a z cosO r ). The simplicity of hj 
comes at the price of introducing an additional gauge 
field = —ihU^ g V Y U g in the transformed kinetic energy 
ftjy2S. The gauge field depends only on the geome- 
try of the local magnetization at the oxide interface. As 
shown below, A^ acts as a q and momentum-dependent 
effective SOI that can be changed electrically because q 
is tunable by a transverse electric field, as shown in FigQ3 

For clarity, we introduce the scaled variables (denoted 
by a bar) r = r/a, q = ag, and k = ak where a is the 
lattice constant and k is the crystal momentum. The 
scaled energy E, and the scaled exchange energy A m 
read 











\ 


r 

/ 3f 






ky/ 


' i 

-A m 



hj 0.0 




FIG. 2: (Color online) Energy bands: E± corresponds to the 
two energy branches, respectively. When Ef < A m , only the 
low energy band E- is involved in the Fermi contour. The 
arrows represent the effective SOI, qk x a y . The strengths of 
the carrier coupling to the local magnetic order is chosen as 
Ef I A m — 1/2 and the spiral wave vector is q — 2tt/7. 



this SOI vanishes. The dependence on k x resembles case 
of a semiconductor 2DEG in a perpendicular magnetic 
field with the Rashba 27 and Dresselhaus^ SOI having 
equal strengths. In this case, when the magnetic-field 
vector potential is taken into account one obtains a reso- 
nant spin Hall conductance; the spin current is carried by 
a charge Hall conductivity 29 . Such a resonance behavior 
is present for a perpendicular spin polarization. In our 
oxide system, however, all averaged values of spin polar- 
ization vanishes due to a zero average magnetization in 
the original spin basis. 

Explicitly diagonalizing the Hamiltonian we obtain 
the eigenenergies 



E = E/e , A m = J/e with e 



2ma 2 

Then the scaled Hamiltonian we find the expression 



H 



(*V S + |^) 2 + {iVyf 



(2) 



(3) 



For a realistic estimates of the parameters of the 2DEG at 
oxides interfaces we choose the lattice constant a = 5A 
To our knowledge, the effective mass m of 2DEG at oxide 
interface is not yet determined. For insulator however, m 
is usually quite large, e.g. for SrTiOa m is ~ 100 times 
larger for GaAs 2 . Here we choose m/m e = 10 with m e 
being the free-electron mass, which sets the unit of energy 
to eo w lbmeV . The Hamiltonian (|3]) we can rewrite in 
the form- 



ic = k 2 x + kl + qk x d y + A m a z 



(4) 



This relation reveal the existence of a SOI that depends 
linearly on q and k, for the collinear spin phase (q — > 0) 



£±(k) = 11 
with the eigenstates 



Am + (qh) 2 



where 



-ikr 



I sin ■ 



ikr 



I sin ^ 

COS -k 



tan 6 = 



qkx 

Am 



COS0 = 



(5) 



(6) 



(7) 



Due to the effective spin-orbit coupling, the Fermi con- 
tours are not parabolic but anisotropic having x and y as 
the symmetry axes, as depicted in Figj2j Although the 
spin states in Eq.© are not independent of k, we still 
have a disappearance of the B erry phase just as the case 
without magnetic field in Ref.[30], which implies that a 
spin current along spin z direction does not exist in the 
absence and presence of an electric field. 
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FIG. 3: The persistent spin current as a function of Ef for 
positive (a) and negative (b) helicities when Fermi level inter- 
sects only the low energy band. The parameters are chosen 
as eo — 15meV, Jo — aeo/2, A m — eo and q — 2tt/7. As 
illustrated in the insets, depending on the direction of an ap- 
plied transverse electric field E z the spiral helicity and hence 
the spin current directions are reversible due to the magneto- 
electric coupling. 



III. PERSISTENT SPIN CURRENT 

The expectation value of spin polarizations per elec- 
tron evaluated using the eigenstates ([6]) are 



sqk x 



VAl + (qk x )2 



(8) 
(9) 



Here s = =b, the double sign corresponds to the two 
branches of the energy dispersion Eq.([5]). In the original 
spin space, (cr|) corresponds spiral spin ordering induced 
by the exchange interaction between the 2DEG and the 
local magnetic moments at the oxide surface. Obviously, 
the y spin-polarization component (a y ) is odd in k x and 
it vanishes upon summing over all occupied states. The 
spin current in the x direction is however generally finite 
when the Fermi level intersects only one of the two bands. 
To prove this we consider the spin current operator, de- 
fined as 



J- 



h 



VjCJi) 



(10) 



where the velocity operators at each k are given by 
v x dH/dp x 2^(2^ + qd y ) and v y dH/dp y 

Considering the symmetry E s (k) = E s (— k) of the 
eigenenergies (Eqj5]), it follows that only (J|) is finite 
and is determined by (in unit of Jo = -p—) 



where f[E s (k)] is the Dirac- Fermi distribution, and k 
runs all occupied states. Introducing the wave- vector pa- 
rameterization k = /c(cos (/?, sin (/?), it can be analytically 
shown that (J|) = when Ef > A m . However, as shown 
in Figj3j the spin current is finite when only the low en- 
ergy band is intersected by the Fermi level. Interestingly, 
the spin current is related to the electron density, rif 
through the Fermi energy Ef. More important, the key 
factor is the odd relationship between the spin current 
and the geometrical spiral structure of the magnetic or- 
dering, being clockwise (q < 0) or anticlockwise (q > 0) 
(Figj3j). This is insofar important, as spin-polarized neu- 
tron scattering experiments 9 on multiferroics evidently 
show that the helicity of the spiral magnetic order is con- 
trollable by a small (~ lkV/cm) transverse electric field, 
as illustrated in Figj3l 

In the absence of an electric field, when the exchange 
interaction is strong enough, i.e. for large A m , the spins 
of the conduction electrons are initially aligned locally 
parallel to n r at each site, (a z ) = 1. Using the Heisenberg 
equation of the electron-spin motion^ we find in the 
linear response regime that the 2DEG develops a uniform 
spin polarization 



(°x) = -q 



eaE x 
2eoA2 



(12) 



when an external electric field is applied along x direc- 
tion. Transforming back demonstrates that the solution 
for (a x ) corresponds to emergence of a spiral spin-density 
wave in the 2DEG rotating in the x — z plane. The di- 
rection of the spin polarization is orthogonal to the oxide 
local magnetic moment. Furthermore, the linear depen- 
dence on q of (a x ) allows for an electric-field control of 
the induced spin helicity. 



IV. HALL CONDUCTIVITY 

The oxide magnetic order is usually not exactly copla- 
nar in the y — z plane but it has a small deviation. 
Here we simulate this non-coplanar modulation with a 
slowly varying spiral order with a spin helicity given by 
(0,/3<2, 0) (/3 <C 1) along y direction. This results in an 
another effective spin-orbit coupling term ~ /3qk y a x with 
a strength (3q. In analogy to the semiconductor case, this 
amounts to the Rashba and the Dresselhaus SOIs having 
different strengths. Therefore, we expect in our case the 
existence of a Hall effect. To see this, we diagonalize the 
resulting total Hamiltonian using the transformation 



qan^&g*.) ) (13) 



where 



<4 y > = 



E 



(2tt) 2 



(2k x a s y +q)f[E s (k)] (11) 



COS( 



, F- k =qJkl + ^kl (14) 
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The Hall effect in the 2DEG is related to the nontrivial 
topology of the resulting eigenstates \k) in the momen- 
tum spac e - 11 ! 19 ! 32 , expressed through the gauge connec- 
tion Afc = —i(k\\/^\k). The off-diagonal Hall conduc- 
tivity is related to Berry's curvature^ f2f = V^; x Aj. 
pointing along the z axis, for which we obtain 



2 cos 2 <p + j3 2 sin 2 ip k dk 

0,1 diverges along the y axis at very small /?, and is sin- 
gular at the origin k = 0. The geometrical Berry phase 
factor 7 S is given by the integral of the curvature over all 
wave vector 




(17) 

The Fermi wave vector k£ is given by the Fermi energy 
Ef = k 2 ± y A 2 ^ + F| . At zero temperature, the off- 
diagonal Hall conductivity a xy for a full band is equal to 
the integral over the Brillouin zone of the component of 
the Berry curvature parallel to z and is thus proportional 
to the Berry phas o 11 ! 19 , i.e. 

<= e ^I n ^=4^ 1 - 1 ^ (i8) 

For a quit small /?, I(/3,kf) — > 0, and a xy = — |^ is 
quantized when the only one of two bands is intersected 



by the Fermi level. Generally, cr xy is not quantized, but 
the transverse conductivity should still be observable. 



V. SUMMARIZING 

A persistent spin current emerges in 2DEG at the in- 
terface of a helimagnet due to the spiral geometry of the 
local magnetic order. The spin current is an odd func- 
tion of the spin helicity and hence electrically controllable 
by a small transverse electric field that reverse the spin 
helicity, making thus a link between spintronics and ox- 
ide electronics. For an in-plane electronic field along the 
spiral we we predict the buildup of carrier spiral spin 
density wave. The spin Berry phase induced by a chiral 
magnetic texture in a Kagome lattice has been discussed 
in Refs.[20|]. Due to a nonzero spin chirality defined as 
the mixed product of three spins on a certain plaque- 
tte, Xijk — Si • (Sj x Sfe), they showed that the Berry 
phase contribution to the Hall conductivity is quantized 
for some values of the band filling. We also calculated 
the Berry curvature and obtained a finite Hall conduc- 
tivity for even a small derivation from the coplanar ox- 
ide helical magnetic order. The transverse conductivity 
is determined by the chirality (/?), the electron density 
and the strength of the exchange interaction (A m ), 
(j xy can thus be quantized, or posseses a nonmonotonic 
behavior upon varying these dependent parameters. 

This research is supported by the DFG (Germany) 
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